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Abstract 

We introduce a new formulation of baryon chiral perturbation theory which improves the ultra- 
violet behavior of propagators and can be interpreted as a smooth cutoff regularization scheme. It 
is equivalent to the standard approach, preserves all symmetries and therefore satisfies the Ward 
identities. Our formulation is equally well defined in the vacuum, one- and few-nucleon sectors of 
the theory. The equations (Bethe-Salpeter, Lippmann-Schwinger, etc.) for the scattering ampli- 
tudes of the few-nucleon sector are free of divergences in the new approach. Unlike the usual cutoff 
regularization, our 'cutoffs' are parameters of the Lagrangian and do not have to be removed. 

PACS numbers: 11.10.Gh,12.39.Fe 
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I. INTRODUCTION 



Weinberg's work in 1979 [l| originated effective field theory (EFT) as one of the most 
important theoretical tools for investigating strong-interaction processes in the low-energy 
regime. The key progress due to Weinberg's approach was the development of a perturbative 
scheme not in terms of a coupling constant, but rather in terms of external momenta and 
the pion mass 1]. In the traditional sense effective field theories are non-renormalizable 
theories. However, as long as one includes all of the infinite number of interactions allowed by 
symmetries, from the point of view of removing divergences there is no difference between the 
so-called non-renormalizable theories and renormalizable theories . Infinities encountered 
in the calculation of loop diagrams are removed by a renormalization of fields and the infinite 
number of free parameters of the most general effective Lagrangian. 

The ideas of Weinberg were further developed and comprehensively applied to the vacuum 
sector of QCD by Gasser and Leutwyler in Refs. Chiral perturbation theory (ChPT) in 

the mesonic sector has been successfully applied to calculations of various physical quantities 
(for a recent review see, e.g., Ref. 0). Processes involving one nucleon in the initial and final 
states were first considered by Gasser, Sainio, and Svarc p. They observed that higher-loop 
diagrams can contribute to terms as low as 0(q 2 ), where q generically denotes a small ex- 
pansion parameter such as, e.g., the pion mass. This problem has widely been interpreted as 
the absence of a systematic power counting in the manifestly Lorentz-invariant formulation 
of baryon chiral perturbation theory (BChPT). As an alternative the heavy-baryon formu- 
lation (HBChPT) was suggested jj], |8| . Most of the calculations in the one-baryon sector 
have been performed in this framework using dimensional regularization in combination with 
the modified minimal subtraction scheme (MS) of ChPT (for an overview see, e.g., Refs. 
El) The advantage of this approach is that it leads to a straightforward power counting. 
Meanwhile it has been realized that, choosing an appropriate renormalization condition, one 
can restore the power countin g in the original manifestly Lorentz-invariant formulation of 

BChPT liEliiElillilliEMlHlSQ- 

A generalization to the few-nucleon sector was sug gest ed in Weinberg's papers on con- 
structing nuclear forces from effective field theory |20l. |21|. For processes involving N > 1 
nucleons, Weinberg proposed applying the power counting to the potential, which is de- 
fined as the sum of all iV-nucleon-irreducible diagrams. The scattering amplitudes are then 
calculated by solving the Lippmann-Schwinger (LS) or Schrodinger equation. 

The application of these ideas has encountered various problems. They originate from 
the renormalization of the LS equation with non-renormalizable potentials (i.e. the itera- 
tion of the potential generates divergent terms with structures which are not included in 
the original potential). A consistent subtractive renormalization requires the inclusion of 
the contributions of an infinite number of counterterms which, in most cases, turns out 
to be technically unfeasible. As a practical solution of the problem one can perform the 
calculations in cutoff EFT. This approach reproduces the results of the subtractively renor- 
malized theory to a given order, provided that the value of the cutoff parameter is suitably 
chosen 




ous applications 



While this approach has been successful in vari- 
the applied cutoff regularization scheme breaks cer- 
tain symmetries of the theory and therefore special care has to be taken. The application 
of cutoff regularization schemes to effective theories has been of interest for a long time 

preserving lattice re gula rization of ChPT in 



the presence of at most a single baryon has been considered in Ref. 43]. Although this reg 
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ularization could, in principle, also be applied in the few-nucleon sector, to the best of our 
knowledge the question of preserving symmetries in calculations of few-nucleon processes 
still remains open. 1 Therefore, the construction of a symmetry-preserving formulation of 
BChPT which renders equations free of divergences is of great interest. 

In this work we use an old idea by Slavnov |44| who introduced chirally invariant terms 
with higher derivatives as a regulator of the non-linear sigma model. We include symmetry- 
preserving higher-derivative terms in the effective Lagrangian of baryon chiral perturbation 
theory which modify the ultraviolet behavior of the pion and baryon propagators. To reg- 
ularize the still remaining infinite number of primitively divergent diagrams 44]| we apply 
dimensional regularization. This ensures that all loop diagrams are regulated. The advan- 
tage of this approach is that it can be applied to individual Feynman diagrams as well as to 
equations of the few-nucleon sector. 

Our work is organized as follows. In Sec. |H] we provide the terms which we add to the 
standard effective Lagrangian. The nucleon mass is calculated within our new approach in 
Sec. 11111 In Sec. II VI we demonstrate that the new scheme satisfies the U(l) Ward identity, 
while in Sec. |V]it is shown that in HBChPT, analogously to the manifestly Lorentz-invariant 
formulation, the existence of a consistent power counting depends on the applied renormal- 
ization condition. Sec. I VII considers an application to simple examples of the nucleon-nucleon 
scattering problem. A summary is given in Sec. IVHI while the Appendix contains the ex- 
pressions for the required loop integrals. 



II. THE MODIFIED EFFECTIVE LAGRANGIAN 

The standard effective Lagrangian consists of the sum of the purely mesonic and the ttN, 
NN, etc. Lagrangians, respectively, 

£ C fT = £tt + CttN + £-NN + • • • • (1) 

The terms in Eq. (H]) are organized in a (chiral) derivative and quark-mass expansion 0, 
11 3 IS H ISM E3- Counting the quark-mass term as 0(q 2 ) 0, 0, the 



mesonic 

Lagrangian contains only even powers, whereas the baryonic Lagrangian involves both even 
and odd powers due to the additional spin degree of freedom. We choose to not show the 
counterterms explicitly. Instead we accompany the Feynman rules with the subtraction rules 
within a fixed renormalization condition. In particular, we use the extended on-mass-shell 
(EOMS) renormalization of Ref. 

The lowest-order mesonic Lagrangian reads 

F 2 

C 2 = — Tr 



DM (IW) f ] + — Tr ( x rf + U X ] ) , (2) 



where U is a unimodular unitary (2 x 2) matrix containing the Goldstone boson fields. The 
covariant derivative is defined as 



DM = d u U - irM + iUl 



ft! 



1 So far, dimensional regularization has only been used for a very restricted number of cases when the 
equations are exactly solvable. 
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where 

r^ = v fJ , + a ll , 1^ = Vf, - a M , x = ^ B (s + ip). 

Here, t> M , a M , s, and p are external vector, axial-vector, scalar, and pseudo-scalar sources, 
respectively. In Eq. (J2J), -F denotes the pion-decay constant in the chiral limit: F v = F[l + 
0(m)] = 92.4 MeV. We work in the isospin-symmetric limit m u = = m, and the lowest- 
order expression for the squared pion mass is M 2 = 2Brh, where B is related to the quark 
condensate (qq)o in the chiral limit 



In the nucleon sector, let 



denote the nucleon field with two four-component Dirac fields, p and n, describing the proton 
and neutron, respectively. The most general ttN Lagrangian is bilinear in $f(x) and \fr(x) 
and involves the quantities u, w M , T^, vj? and X ± (and their derivatives), which are defined 

as 

u 2 = U, u ll = iv)DJJu\ r M = - [v)d^u + ud^v) - i^r^u + ul^v))] , 

X± = u ] xu ] ± ux ] u. 
In terms of these building blocks the lowest-order Lagrangian reads 0, 0] 

= * (il»D» -m+ l -°9A 1,1^ % (3) 

where = (d^ + T fl — iv^)^/ denotes the covariant derivative. (In the definition of the 
covariant derivative we follow Ref. jHj, where only contains traceless external fields and 

( \ - ° 

the coupling to the isosinglet vector field vf is considered separately.) In Eq. (jSJ), m and 9a 
refer to the chiral limit of the physical nucleon mass and the axial- vector coupling constant, 
respectively. 

Below we will calculate the nucleon self-energy to third order. For that purpose, we will 
need one of the seven structures of the Lagrangian at 0(q 2 ) 0, H?]], 

C®= Cl Tr( X+ )W + .... (4) 

(3) 

The Lagrangian C n ^ does not contribute to the nucleon mass at the given order. 

To improve the ultraviolet behavior of the propagators generated by the Lagrangian of 
Eq. (^) we introduce additional terms into the Lagrangian which modify the propagators of 
the pion and the nucleon. In particular, we consider the modified pion propagator 

A - ^ = p 2_ M 2 + , + II A 2 _ 2_ , + ( 5 ) 

j=l 

and the modified nucleon propagator 

1 N * A 2 
(„) = 1 TT (a) 

NyP> (tf-m + i0+) A J AL + m 2 - p 2 - t0+ ' { ' 
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Here, A ni and are (independent) parameters. 2 For simplicity we use the standard 
prescription for dealing with poles in the A-dependent factors of the modified propagators. 
For sufficiently large values of the parameters A any other prescription leads to the same 
results for low-energy physical quantities. The propagators above can be generated by 
a Lagrangian which, in addition to the standard BChPT Lagrangian of Eq. JI}, contains 
additional symmetry-preserving terms. These terms vanish in the limit A^ — > oo, A^- — > oo. 

The choice of the additional terms of the Lagrangian is not unique. Furthermore these 
terms not only generate the above propagators, but also result in additional interaction 
terms. Our choice is motivated by the simplicity of calculations. For the pion sector we 
choose 



F 2 



n=l 



4 4 



(D 2 ) n UU^ - U [(D 2 ) n U] D 2 UW - U (D 2 Uy - X U ] + U X 



where D 2 U = D a D a U and X n are functions of A^. For example, in order to generate the 
modified propagator 



Aj(p) 



we need to take N n = 3 and 



n 



A 2 - 



p 2 — M 2 + i0+ 11 A 2 , - p 2 - i0+ 

.7 = 1 nj 



(7) 



x 2 
x 3 



1 



1 



1 



A 2 A 2 A 2 ' 

iV 7Tl IV tt2 IV tt3 
A2 i A 2 _| — A2 



A 2 A 2 A 2 

1 



A 2 A 2 A 2 ' 

For the additional terms of the Lagrangian of the nucleon sector we choose 

A'vp 



(8) 



Civ = Ey - m) (D 2 + m 2 )" * + h.c] , 

n=l 

where Y n are functions of A^. For example, for the modified nucleon propagator 



(9) 



A 2 



- m + i0+) (A| + m 2 - p 2 - i0+) 



(10) 



we have Nq, = 1 and Y\ = 1/A|. 

Depending on the order of the performed calculations we choose the modified propagators, 
i.e. fix N n and Ny, such that all loop diagrams (except some of the primitively divergent 
diagrams) contributing to the given order converge. To obtain the modified Lagrangian for 
lower order calculations one needs to take A^ — > oo, A$j — > oo for some of the parameters 
A in the modified Lagrangian used in higher-order calculations. 



In the following we let A collectively represent the and A*j. 
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Analogously to the nonlinear sigma model 44]| the additional terms do not render all loop 
diagrams finite. There still remain an infinite number of primitively divergent diagrams in 
the mesonic sector as well as divergences in diagrams with fermion loops. These diagrams 
can be regularized in a symmetry-preserving way by introducing additional auxiliary fields 
analogously to the case of Yang-Mills theory j53|. However, in practical calculations such a 
technique is rather difficult to apply. Instead it is possible (and much more convenient) to 
use standard dimensional regularization. This is due to the fact that the remaining divergent 
diagrams contribute either in physical quantities of the vacuum (purely mesonic) and the 
one-nucleon sectors, or they appear as sub-diagrams in the potentials of the few-nucleon 
sector. In both cases the calculations are perturbative, i.e., to any given order in the chiral 
expansion one needs to calculate a finite number of diagrams. Therefore, divergences which 
show up as the l/(n — 4) poles (where n denotes the number of space-time dimensions) can 
be explicitly subtracted (i.e. absorbed in the redefinition of the parameters of the effective 
Lagrangian). 

To summarize, our scheme consists of adding symmetry-preserving additional terms in 
the standard effective Lagrangian and applying dimensional regularization to the resulting 
effective theory. All symmetries are preserved in the regularized theory, i.e., regularized 
quantities satisfy all relevant Ward identities. We expand the regularized diagrams in powers 
of n — 4 and subtract l/(n — 4) pole-terms observing that there is a finite number of them to 
any given (finite) order in the chiral expansion of physical quantities in the vacuum and one 
nucleon sectors, and the potentials in the few-nucleon sector. No further divergences occur 
(for finite parameters A) neither in the vacuum and one-nucleon sector nor in the equations 
of the few-nucleon sector. Therefore, for the equations of the few-nucleon sector we can take 
n = 4. 

Using a field transformation the additional higher- derivative terms which we introduced in 
the effective Lagrangian can be reexpressed in a canonical form, i.e. a form with a minimal 
number of independent terms [3, 55[. This clearly shows that any A dependence of the 



physical quantities can systematically be absorbed in the redefinition of the parameters of 
the standard canonical effective Lagrangian. 



III. NUCLEON SELF-ENERGY 

As an example of the application of our approach we calculate the nucleon self-energy to 
order O (q 3 ) in this section. For this calculation it is sufficient to take N n = = 1. We 
parametrize the complete nucleon propagator as 

Sn ^ = J - m + i0+) [1 - {p 2 - m 2 ) /A|] - Sp N ' ^ 



where m is the nucleon pole mass in the chiral limit and the nucleon self-energy — iE( 
represents the sum of all one-particle-irreducible perturbative contributions to the two-point 
function. The physical nucleon mass is defined through the pole of the full propagator at 
i = fn N , 



{mN _ m) I i _ m N m ) _ S(m7V ) = o. (12) 



A 



2 

(2) 



At 0(q ), the self-energy receives contact contributions from as well as the one-loop 
contribution of Fig. [TJ 

+ Si 00 p, (13) 



J contact 
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FIG. 1: One-loop contribution to the nucleon self-energy at 0(q i ). 



where S conta ct = — 4ciM 2 . Applying Feynman rules we obtain for the one-loop contribution 
Sioop = - 3 ^ A * {(*+™) /(1011) + M 2 (^+m) 7(1111) 



4F 2 



+ (p 2 -m 2 )p l J (p) (llll)-^ J (p) (1110)} 



where 



with 



{i (abed), p»I {p) (abed)} 



d 4 k {l,k"} 
(2vr) 4 A a B b C c D d '' 



(14) 
(15) 



A = k 2 - A 2 + z0 + , 

B = k 2 — M 2 + i0 + , 

C = (p + A;) 2 - m 2 - A% + z0 + , 

D = (p + k) 2 — m 2 + i0 + . 

To further simplify the calculation we take = A#. The parameter n of dimensional 
regularization has been put n = 4, as the diagram is finite for finite A. We perform the 
renormalization by applying the extended on mass- shell (EOMS) scheme of Ref. 0. First 
we substitute the expressions for the loop integrals from the Appendix and expand Eq. (JTljl 
in a power series in A (around A = oo). We then subtract all positive powers of A and 
ln(A/m). 3 Next we expand the remaining expression in powers of small quantities, i.e., M, 
p 2 — m 2 and f 1 —m and subtract all terms of zeroth, first and second order in this expansion, 
so that the renormalized expression is indeed of order q 3 as mandated by the power counting. 
The resulting expression for the subtraction terms reads 



o2 



-isub 



3 9 A A 2 (4m + hp 1 ) 
256tt 2 F 2 



+ 



o2 

9 A 



2567r 2 mF 2 



„2\2 



-12m 4 + 3m M + 8m V ~ Hp ) 



A 



•10m 3 $f + mp 2 ^ + 12m (2m 3 + 4mM 2 + 3m 2 ^ - p 2 ^) In ( — 



in 



o2 

9 A m 



2560tt 2 F 2 A 2 



-48m 4 - lOQo 2 ) 2 + 6m 2 (15M 2 + 41p 2 ) + 63m 3 jrf + UOmp 2 ^ 

A 



240m 2 (m 2 + 2M 2 + p 2 + 2m#S) In 



m 



(16) 



3 Note that, since our scheme respects all symmetries of the theory, the Ward identities are satisfied sepa- 
rately in each order of the expansion in powers of A. 
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Subtracting Eq. (|16p from Eq. (|14p and taking f 1 = tun, we obtain for the renormalized 
on-mass-shell self-energy to order q 3 



^ \p=m N 



3g A M 3 



32vrF 2 

Using Eq. , the nucleon mass to order q 3 follows as 



rriN = m — Ac\M 



3 9 A 
32ttF 2 



(17) 



which agrees with the standard BChPT result [3 E, H ■ 



IV. ELECTROMAGNETIC WARD IDENTITY 

To demonstrate that the new formulation indeed respects the symmetries of the theory, we 
analyze the electromagnetic Ward identity for the nucleon which, in units of the elementary 
charge, reads 

( Pf - Pl )^% (p /lft ) = [S^ (p f ) - S N ' (j^)] . (18) 

Here, 

(pf, Pi) = T% (p f , pi) + A% (p f , Pi ) (19) 

is the one-particle-irreducible three-point function (\&j7 M ^ r ) with the electromagnetic 
current operator in units of the elementary charge. F^ (Pf,Pi) corresponds to the tree- 
order contribution and A^(pf,Pi) consists of loop corrections. In order to determine T^, 
we consider the coupling to an external electromagnetic field A^ and insert for the external 
fields in Eq. © 

For the purpose of this section it is sufficient to take N n = = 1. From our modified 
Lagrangian we obtain 

(Pf,Pi) = ^ 7* ~ ^ [V (Pf + Pi - 2^ 2 ) + (Pf +PiT (h + h ~ 2m)] , 



(20) 



2 2 2A 2 

where r^ (pf,Pi) and the free propagator of Eq. (|1(J|) satisfy the relation 

(Pf ~Pi)^ m fa,*) = IsST 1 (Pf) - S^ 1 (Pi)) . (21) 



Of course, this result is not surprising, because the coupling to an external electromagnetic 
field in the Lagrangian of Eq. (jOJ) proceeds via covariant derivatives which essentially amount 
to a minimal coupling. At tree level this automatically results in contributions satisfying 
the Ward identity. 4 



4 In the context of EFT the use of minimal-substitution terms alone is not sufficient to generate a consistent 
framework, because the most general effective Lagrangian also contains terms involving field-strength 
tensors such as, e.g., the l§ and Iq terms of £4 Jj. In general, the presence of these terms is also necessary 
for the purposes of renormalization (see Ref. |58| for a critical discussion of this issue) . 
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(a) 



-►- + 



Pi Pf 
(b) 



— + 



X 



Pi 



Pf Pi 



Pi 



(c) 



(d) 



FIG. 2: One-loop contributions to the electromagnetic vertex. 



For the one-loop corrections to the nucleon self-energy (of Fig. [I} and the vertex (diagrams 
of Fig. |2J) we obtain 



o2 



(22) 



A^ a (p f ,Pi) 



t A . f d n k 
AF 2% J (27r) n 
xr^ ( Pf + k,p t + k) S± (ft + fc) ^75 r a AA(Ar), 



o2 



cTA; 



A£ 6 (p /)Pi ) = 2r 3 ^ y -^^75 fift( P/ + A0 7 *76Aj(*) l 



o2 



Kc(Pf,Pi) = 2r 3 -^i 



2ttJ 
d n k 

AF 2 ' J J2^Y 
3 A . /" cTfc 



xr; (p, - Pi - k, -k) 1{ 75 A£(fc + p, - p/)A£(A;) ; 



where 



K (p',p) = (p'+pr(i 



p' 2 +p 2 - M' 

a! 



(23) 
(24) 

(25) 

(26) 
(27) 



\ 7T / 

is the leading tree-order contribution in T^, which is related to the one-particle-irreducible 
three-point function (iTjJ' fl 7i i ) by the relation 

I% i tf,p)=iewl»tf,p). (28) 

To check the Ward identity for the above loop diagrams, we multiply = A^ a + A^ b + 
A^ c + A^ d with (j)f — Pi) , use Eq. (|2~Tj) and the Ward identity for pions (at leading tree 
order) 

(p'-p),Ko(p',p) = Af V) " Aj-'Cp), (29) 
and obtain after a straightforward calculation 

(Pf ~ Pi)^K (Pf,Pi) = s ioop(^J - Si oop (^ f ) , (30) 



which verifies the Ward identity of Eq. ()18|) . 



9 



V. NUCLEON SELF-ENERGY DIAGRAM IN HBCHPT 



It is common practice to assume the existence of a consistent power counting in HBChPT 
without specifying the renormalization scheme used. In HBChPT, as in any quantum field 
theory, one has the freedom to choose a renormalization condition. Dimensional regular- 
ization in combination with the MS scheme, which is commonly used in HBChPT, is only 
one among an infinite number of possibilities. In this section we apply our higher-derivative 
formulation to the nucleon self-energy diagram of Fig. [I\m order to show that, analogously 
to the manifestly Lorentz-invariant formulation, the existence of a consistent power counting 
in HBChPT depends on the choice of the renormalization condition. 

Using the pion propagator of Eq. (jSJ) for N n = 2 and A„-i = A w2 = A, we obtain (see, e.g., 
Section 5.5.9 and Appendix C.2 of Ref. [5J] for a detailed calculation in HBChPT) 

° 2 4 

S&CP) = 3 %r S^S: J&(121; u), (31) 

where uo = (p ■ v — m) and 

J^r(abc:uj) — % f - — — ; . (32) 

nNK J (27r) 4 [£;2_M 2 + 20+] a [P-A 2 + 20+]V k + cu + tO+Y 

One can parameterize J^-(121;u;) as 

J^ N {\21;u3) = c x gr + c^v v . (33) 

Since S v ■ v = 0, c 2 does not contribute to the self-energy. For c\ we find 

ci = \ [{M 2 - u 2 ) 7^(121; u) + J^(021; u) + to J wJV (120; u)] , (34) 

where 

/d 4 k 1 
7- r . (35) 
(2vr) 4 [jfe2 _ M 2 + iQ +] a [jfc2 _ A 2 + iQ +] b [|, • k + UO + i0+] C 

Standard power counting of HBChPT assigns the order O (q 3 ) to the diagram of Fig. 
Calculating the loop integrals of Eq. (see Appendix |Bj), we obtain 

2 

E £op (P) = [- A3 + 3 ^ A 2 ] + O(A). (36) 

Both terms inside the square brackets (as well as the term proportional to A which, for the 
sake of brevity, we have not displayed) violate the power counting. They are analytic in 
momenta and can be absorbed in the renormalization of the nucleon mass and the nucleon 
field. Note that the corresponding mass counterterm SmN v N v , cancelling the A 3 term in 
Eq. ([36)1 . is equal to zero in the standard formulation of HBChPT with dimensional regu- 
larization and is, therefore, usually not indicated in the effective Lagrangian of HBChPT. 5 
Choosing the renormalization scheme appropriately one can subtract all terms in Eq. (|36|) 
which violate the power counting so that the renormalized diagram is of order O (g 3 ) . 



5 This is analogous to the case of the pion t adp ole self-energy in cutoff regularization, where one needs a 
counterterm of order p 2 for the pion mass | |52| . 
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+ 




FIG. 3: Graphical illustration of the equation for the NN scattering amplitude. 

VI. NN SECTOR 

A. Contact interaction 

In this section we consider a demonstrating example of the application of our approach 
to the NN problem in the manifestly Lorentz-invariant formulation of baryon chiral pertur- 
bation theory. 

Let us consider the simplest possible (\i/\I/) 2 contact interaction term. The corresponding 
equation for the lowest-order amplitude can be solved analytically, therefore, one can apply 
standard dimensional regularization to this problem. The interaction term in the Lagrangian 
reads 

C NNNN = C M M. (37) 

In the following, we will consider the scattering amplitude in the center-of-mass frame. 
Let P = pi + P2 denote the total four-momentum of the scattered nucleons, where pi = 
( a/ m 2 + p 2 ,p), P2 = {\/m 2 + p 2 , —p) with p = \p \ and P 2 = Am 2 + Ap 2 . 

The interaction Lagrangian of Eq. (J37)) generates the following NN vertex (two-nucleon- 
irreducible contribution in the NN scattering amplitude) 

iV\c^v = 2iC (SxuSa^ - S XfJ ,S av ) . (38) 

According to Weinberg's approach, to find the corresponding lowest-order NN scattering 
amplitude we need to solve the equation 6 

/d 4 k 
yj Vao-,0!7 Gaj,ps(Pjk) Tps,fiv(P), (39) 

schematically shown in Fig. El where 



Gfxy^siP, k) 



[k 2 -m 2 + i0+] [(P - kf -m 2 + z0+] 



A 4 

G aj ,p S (P,k) = - ^ - . q+] ^ p *_ — 2 _ m2 _ A| + ^ G ajtPS (P,k) 

are the (two-nucleon) propagators to be used in standard dimensional regularization and 
higher-derivative formulation, respectively. Integrating Eq. ()39j) over k, we obtain 

T\ a ^ u (P) = V X c,iw + iV X a,ay Gay,/3S(P) T^ V {P) , (40) 



It is understood that T\ a ^ v needs to be multiplied with the corresponding Dirac spinors. 
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where 



z)orf,j38 (P) = i^rnSap (rf s P(i + mdys) Inn 

+ [lap (ljs p v + m5 l5 ) - mSaptfg] 4 P ^P M - i^ p i u lS lNN,nu \, (41) 



with 



Inn, P 1 NN , 1 nn $ -tj { ^ )n _ ^ + _ 2 _ 



(42) 



(2tt)" [A; 2 - m 2 + zO+] [(P - fc) - m 2 + iO+] 
in standard dimensional regularization and 

I NN ' NN ' NN i ' J (2vr) 4 [k 2 - m 2 + .0+] [(P - k) 2 - m 2 + *0+] 

A 4 

x F^ 1 2 T" (43) 

[Jfc2 _ m 2 _ A 2 + i0 +] [(p _ A;) 2 - m 2 - A| + z0+] 

in higher- derivative formulation, respectively. 

We renormalize Eq. (|4()|) by subtracting the contributions of loop integrals at P 2 = 4m 2 . 
Next we expand the subtracted loop integrals (I R = I — J|p2 =4m 2) in p and retain terms to 
order O (p). The resulting equation reads 

T^ V {P) = Vx^u + iV Xa , a , G^8 7>V(P), (44) 

where 

P PM pv \ 

+ f [-fop fe^ + ™M " ~ -faffa —f- j • (45) 



Comparing Eq. (|44|) for standard dimensional regularization and higher-derivative formu- 
lation (remembering that terms of order 0(p 2 ) have been neglected) we obtain identical 
results. 7 

Note that, in contrast to the example considered above, it is not clear how to apply 
standard dimensional regularization to equations involving potentials derived from BChPT. 
The advantage of our higher- derivative formulation is that it is also applicable in these cases 
while preserving the symmetries of the theory. 



B. Inclusion of one-pion exchange potential 

Below we consider some conceptual issues of renormalization and the A dependence of the 
leading-order NN scattering amplitude in a non-relativistic formalism of BChPT. This am- 
plitude is obtained by solving the Lippmann-Schwinger equation with a contact interaction 
plus the one-pion exchange potential. 8 



7 Although Eq. can be solved exactly, it is beyond the scope of this paper to perform this straightforward 
but rather cumbersome calculation. 

8 A detailed discussion of the heavy-baryon reduction of our new Lagrangian (including a numerical analysis 
in the few-body sector) will be given in a forthcoming publication. 
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Using an appropriate field redefinition and the standard heavy-baryon reduction with 
v = (1, 0, 0, 0) we obtain the leading-order NN potential (for the choice = 1, N n = 0) 

V(p,p)=C s + C T a 1 .a 2 -^—j (n • r 2 ) ( ^_. )2 + M , , (46) 

where C5 and CV are the coupling constants of the four-nucleon contact interaction La- 
grangian at leading order. The scattering amplitude satisfies the equation 

T (p',p) = V (p',p) +m f V (p', k) A2 ^ T (k,p) , (47) 

J {2ny V M A 2 + kA (mE - k 2 + i0+) V J 

where E = p 2 /m is the energy of two nucleons in the center-of-mass system. 

We fix the free parameters C,s and CV as functions of A by demanding that the solution 
of equation (jT7|) reproduces two physical quantities for fixed kinematics. Following Ref. [2l[ 
we choose the renormalization points of the order of small external momenta. This exactly 
corresponds to the following general renormalization program (see, e.g., [E^]). First one 
calculates the quantities of physical interest in terms of bare parameters in the regularized 
theory. Once a sufficient number of physical quantities are determined as functions of bare 
parameters, one inverts the results and expresses the bare quantities in terms of physical 
quantities. These expressions are then used to eliminate the bare parameters in all other 
quantities of physical interest. This procedure preserves all symmetries provided that the 
applied regularization scheme respects them. If the considered theory is renormalizable in 
the standard sense the above procedure removes all divergences. 

For definiteness, let us take the zero-kinematics as renormalization points. 9 This corre- 
sponds to the subtraction of loop diagrams at zero-kinematics. In the above case, expressing 
Cs and Ct from two physical quantities 10 and substituting them into other quantities we 
only eliminate some of the terms that diverge in the limit A — > 00. This is due to the non- 
renormalizability of BChPT in the traditional sense. However, note that A is a parameter 
of the Lagrangian and we do not have to take it to infinity. The remaining A dependence of 
the amplitude is of higher order in the small-parameter expansion (pion mass, small external 
momenta). As the potential of Eq. (J46)) is non-renormalizable in the traditional sense, the 
perturbative expansion of the renormalized amplitude contains negative as well as positive 
powers of A (and/or positive powers of In A). These contributions contain terms of the form 

q i q i A j 
~ — as well as ~ —. , . , with i > 0, j > 0, (48) 

where q denotes small external momenta or the pion mass and Q stands for AirE and/or the 
large scale parameter hidden in renormalized contact interaction constants. To keep these 
formally higher-order contributions indeed suppressed numerically, one should take A ~ Q. 
The existence of such an optimal value of the parameter A depends on the validity of the 
assumption of Weinberg's approach that the renormalized coupling constants are natural 



9 This would not be a good choice if we took the A — > 00 limit in the end. We would be faced with the 
problem of very poor convergence jjof . 
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We could take as "quantities of physical interest" the scattering lengths of the 1 S'o and 3 Si NN scattering. 
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for renormalization points of the order of or less than small external momenta. The validity 
of this assumption has to be checked at each order of calculations. While one cannot take 
the existence of the optimal value of A for granted, the reasonable success of cutoff EFT 
suggests that it should exist. A detailed analysis of this issue in our symmetry-preserving 
approach is in progress. 

The complete A dependence of physical quantities can be absorbed in the redefinition of 
couplings constants of the canonical Lagrangian. 11 For the above example this means that 
the contributions of the form of Eq. ()48|) can be absorbed in the redefinition of higher-order 
coupling constants. We could take any value for the parameter A provided that the compen- 
sating contributions of higher-order terms (an infinite number of them) are also included, 
but that does not seem to be feasible. Note that A is not a cutoff-regularization parameter 
and does not need to be taken to infinity. The above specified optimal choice of the free 
parameter of the Lagrangian, A, ensures that, to the accuracy of the given calculations, phys- 
ical quantities do not depend on higher-order terms which we introduced in the Lagrangian 
(A independence). 

As one cannot solve equations exactly, one carries out the above renormalization program 
numerically by fixing coupling constants as functions of A so that the given particular phys- 
ical quantities at the renormalization points (i.e. the fixed A-independent values of them) 
are reproduced. The reliability of this numerical renormalization procedure in comparison 
with the explicit analytic renormalization depends only on the accuracy of the numerical 
approximation, i.e. the two approaches are conceptually equivalent. 

In the approach suggested in this work, Ward identities are satisfied order by order 
in the loop expansion as well as in the chiral expansion of physical quantities 61] . In 
the few-nucleon sector the physical quantities, at any finite order in the chiral expansion, 
contain an infinite number of terms in the loop expansion. On the other hand, to any 
specified order q n , for a processes involving A nucleons, there is a finite number of A- 
nucleon irreducible diagrams. The sum of these diagrams is defined as the effective potential. 
An infinite number of diagrams contributing in the scattering amplitude (at given order 
q n ) is summed up by solving the corresponding equations with given n-th order effective 
potential. I.e. substituting the g n -order potential in the Lippmann-Schwinger equation and 
performing the renormalization properly (as specified above) corresponds to the summation 
of all renormalized diagrams up to order q n . The solution of the Lippmann-Schwinger 
equation also contains some, but not all, of the higher-order contributions, and the result is 
reliable only to order q n , the error being of order 0(q n+1 ). As the Ward identities are satisfied 
order by order in the chiral expansion, and the Lippmann-Schwinger equation resums all 
contributions to order q n , the contributions in physical quantities which can violate the 
identities are of order ~ 0(q n+1 ), i.e. beyond the accuracy of the given calculations. 



11 The original coupling constants are written as c, = c[ + Sci, where the c[ are redefined coupling 
constants and the loop expansion of the 5ci part exactly cancels the corresponding A-dependent parts of 
loop diagrams. The c[ are independent of momenta, i.e. local interaction terms of the effective Lagrangian 
remain local. 
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VII. SUMMARY 



We have discussed a new formulation of BChPT, which preserves all symmetries of the 
theory. The main idea is to use some of the structures of the most general effective La- 
grangian to improve the ultraviolet behavior of propagators. The coefficients of these terms 
depend on parameters (with dimension of mass) which serve as smooth cutoffs of the theory. 
For practical applications it is convenient to choose these parameters to be equal. 

We have explicitly applied our new approach to a calculation of the nucleon mass to order 
O (g 3 ). We have also explicitly verified that the electromagnetic Ward identities are satisfied 
by (strong) one-loop-order corrections. The application of this scheme to the one-nucleon 
sector of HBChPT demonstrates that the existence of a consistent power counting scheme 
in HBChPT actually depends on the applied renormalization scheme. 

The considerable advantage of the new formulation in comparison with standard dimen- 
sional regularization is that, while preserving all symmetries of the effective theory, it leads to 
equations in the few nucleon (NN, NNN, etc.) sector which are free of divergences. We have 
explicitly considered examples of the contact interaction and one-pion exchange potentials 
in the NN scattering problem and have discussed issues of renormalization and consistency. 
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APPENDIX A: NUCLEON SELF-ENERGY 

The explicit expressions for the loop integrals of Eq. (|15|) contributing in the calculation 
of the nucleon self-energy up to and including order O (-h) read 
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APPENDIX B: HBCHPT 

The expansions of the considered heavy-baryon integrals of Eq. (|33j) around A = 00 are 
riven by 
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APPENDIX C: NN SECTOR 

1. Standard dimensional regularization 

The explicit expression for the loop integral Inn in dimensional regularization is given 

Inn = 2A + I NN , (CI) 
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2. Higher-derivative formulation 

In the following, the explicit expressions for the loop integrals in higher-derivative for- 
mulation up to and including order O are given. The scalar integral reads 
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and, after subtraction, one obtains 
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